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At the Karlsruhe pulsed 3.75 MV Van de Graaff accelerator the ^*'S(n,7)^''S(5.05 min) cross 
section was measured by the fast cyclic activation technique via the 3.103 MeV 7-ray line of the 
1/^ ' ^^S-decay. Samples of elemental sulfur enriched in '^^S by 5.933 % were irradiated between two gold 

, foils which served as capture standards. The capture cross section was measured at the neutron 

0^ ' energies 25, 151, 176, and 218 keV, respectively. The ^^S(n,7)^^S-cross section in the thermonuclear 

and thermal energy range has been calculated using the direct-capture (DC) model combined with 
the folding procedure used for the determination of the potentials. The non-resonant experimental 
data for this reaction can be reproduced excellently using this method. The input parameters of the 
DC-calculation (masses, Q-values, nuclear density distributions, spectroscopic factors, spin-parity 
1^-^ ■ assignments and excitation energies of the low-lying states of the residual nucleus) have been taken 

' from the available experimental data. 
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^ . PACS numbers: 25.40.Lw, 24.50.+g, 25.40.Dn 

> 

O ■ I. INTRODUCTION 

o . 

jy-^ ' In the last years the importance of the direct-reaction (DI) mechanism in nucleosynthesis has been realized. The 
DI dominates over the compound-nucleus (CN) reaction mechanism if there exist no CN-levels near the threshold 
that can be excited in the reaction. For instance, this can be the case in reactions involving light nuclei, in the big- 
, bang scenario as well as in in stellar hydrogen and helium burning. Direct capture (DC) can also be of importance 
' I in proton capture by proton-rich target nuclei in the rp-process occurring in novae or X-ray bursts and in 

O neutron capture by neutron-rich nuclei [Q- Jill in the inhomogenous big bang as well as in the s-process taking place 
2 in helium burning of red giants. Furthermore, DC can dominate for neutron capture by light ^ and heavy fl^j- 
^ target nuclei far from stability in the a-process and r-process occurring in supernovae of type II, respectively. 

^ , The reaction rates of the neutron-rich S-isotopes are of interest in the nucleosynthesis of nuclei in the s-process 
' i'^ the S-Cl-Ar-Ca region |jl^ , |jl^ , inhomogenous big-bang scenario , [|l^ , and in the a-rich freeze out of the 
i^jj , neutron-rich hot neutrino bubble in supernovae of type II p^ - [pl| . 

For the rare isotope '^^S a significant abundance contribution is expected from the s-process nucleosynthesis. For 
quantitative analyses the size of the destruction rate, i.e. the neutron capture rate, is of fundamental importance 
to estimate the magnitude of the "^^S abundance formed by the weak and main s-process components. Using the 
statistical model the ^^S capture cross section has been estimated to be SOO/ubarn at 30keV by Woosley et al. [ p2| . 
The present measurement, the first experimental investigation of this cross section, applies the fast cyclic activation 
technique developed at the Karlsruhe 3.75 MV Van de Graaff accelerator. 

We investigate the capture reaction •^^S(n,7)'^^S from thermal (25.3 mcV) to thermonuclear (25-218keV) projectile 
energies and compare the calculated cross sections in the DC-model with the experimental data. In Sections || and III 



the activation technique is described, and the experimental results are given. In Section IV we introduce the methods 
used for calculating the direct-capture (DC) cross sections. In Section ^ the experimental and theoretical results for 
the neutron-capture cross section of '^^S are compared and discussed. 



II. FAST CYCLIC ACTIVATION TECHNIQUE 

The measurements have been carried out at the Karlsruhe pulsed 3.75 Van de Graaff accelerator. A common 
activation measurement is subdivided into two parts: (1) the irradiation of the sample, (2) the counting of the 
induced activity p3] . The cyclic activation method is the repetition of the irradiation and activity counting procedure 
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many times to gain statistics. Especially for nuclei with half lives of only minutes or seconds, in our case '^^S with a 
half life of 5.05 minutes, a large number of irradiation and counting cycles is needed. The time constants for each cycle 
which are chosen shorter than the fluctuations of the neutron beam and comparable or shorter than the decay rate A 
of the measured isotope are, the irradiation time, tb, the counting time, tc, the waiting time, (the time to switch 
from the irradiation to the counting phase), and the total time, T—tb+t^+tc+f^ (t'^ the time to switch from the 
counting to the irradiation phase). In the actual '^^S measurements the runs were partly carried out with tb=19.58s, 
tc=19.24s, T=40s, partly with tb=49.58s, ic=49.24s, T=100s. The waiting time is in both cases iw=0.42s. 

The accumulated number of counts from a total of n cycles, C ~ X]"=i where C;, the counts after the i-th cycle, 
are calculated for a chosen irradiation time, i;,, which is short enough compared with the fluctuations of the neutron 
flux with H 

C = t^K^j^\-\\ - exp(~Atc)] exp(-AU)^^^^^^^^^^a[l - A exp(-AT)] (i) 

i exp^ Ai ) 



with 

Er=i*^exp[-(n-i)AT 



The following additional quantities have been defined: e^: Ge-efficiency, K^: 7-ray absorption, /-,,: 7-ray intensity 
per decay, N: the number of target nuclei, cr: the capture cross section, $^1 the neutron flux in the i-th cycle. The 
quantity /b is calculated from the registered flux history of a ^Li glass monitor. 

The efficiency determination of the 35% HPGe-detector (2keV resolution at 1.332 MeV)has been reported else- 
where 0] . The 7-ray absorption was calculated using tables published by Storm and Israel and Veigele |Q . The 
half-lives and the 7-ray intensities per decay of ^^S and ^^^Au are given in Table |. 

The activities of nuclides with half lives of several hours to days, i.e. the activity of ^^^Au, can be also counted after 
the end of the cyclic activation consisting of n cycles using 

n 

C„ = e^if^AA-i[l - exp(-ATM)] exp(-ATw)[l - exp(-Aib)]iVfT/b ^ . (2) 

i=l 

Here Tm is the measuring time of the Ge-detector and Tw the time elapsed between the end of cyclic activation and 
begin of the new data aquisition. 

The Eqs. |l| and 1^, respectively, contain the unknown quantities a and the total neutron flux ^Y^=\ Therefore, 
cross section ratios can be formed for different isotopes exposed to the same total neutron flux. This is the basis for the 
determination of the "^^S capture cross section relative to the well-known standard ^^^Au capture cross section [ p7| . 
As the "^^S sample to be investigated is characterized by a finite thickness it is necessary to sandwich the sample by 
two comparatively thin gold foils for the determination of the effective neutron flux at sample position. The activities 
of these gold foils were counted also individually after termination of the cyclic activation. The effective count rate 
of gold was obtained from these individual rates as well as from the accumulated gold count rate during the cyclic 
activation run. Therefore, the effective neutron flux at sample position was determined in two ways by way of the gold 
activation according to the Eqs. |l| and ^ Using Eq. |l| has the advantage that saturation effects in the gold activity 
for irradiations over several days are avoided ||23| . 



III. EXPERIMENTAL ARRANGEMENTS AND RESULTS 

In Fig. a scheme of the experimental setup is shown. The kinematically collimated neutron beam is generated 
with the ^Li(p,n) reaction near the reaction threshold (1.912 MeV proton energy) with thick ^Li-targets (30 /im) and 
corresponds to a Maxwellian neutron spectrum with a thermal neutron energy of kT = 25 keV [^ p?! . The neutron 
spectra at the neutron energies 151, 176, and 218keVkeV were generated using thin Li-targets (2.5 fj,m). The required 
proton energy conditions and the neutron spectra integrated over the solid angle of the sample were determined in 
time-of-flight (TOE) measurements before the actual activation runs using the accelerator in pulsed mode (Fig. ||). 
These neutron spectra are in good agreement with corresponding Monte Carlo calculations |2^ and |2^ . The proton 
beam was wobbled initiated by magnetic deflection to cover the area of the Li-target. The beam profile formed was 
studied on a quartz target. To switch back and forth between sample irradiation and activity counting a fast sample 
changer operating with compressed air was used. Close to the beam line where the neutrons have been generated the 
Ge-detector for activity counting, well shielded by lead and Li-loaded paraffin, has been installed (Fig. ||). During 
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the irradiation phase the analog to digital converter was gated to prevent data acquisition. The relative neutron 
flux was recorded continuously with a ®Li glass detector. During the activity counting phase neutron generation was 
interrupted by a beam stop for the proton beam. This is essential to reduce all prompt accelerator dependent 7-rays. 
The beam stop was installed in the beam line at the accelerator hall so that in the activity counting periods the 
experimental hall was free of prompt background radiation. In Table | the sample characteristics of '^^S and ^^^Au 
and the decay properties of the product nuclei are listed. Table || gives a survey of the sample weights and the 
measured '^^S capture cross sections. The sulfur was pressed to self supporting tablets of 6 mm diameter, the gold 
foils on back and front side of the sulfur sample with regard to the impinging neutrons had the same dimensions. 
At 25keV neutron energy measurements were carried out with sulfur sample masses between 20 and 100 mg. No 
significant effect from multiple neutron scattering was observed. A sample of the sulfur powder was also heated to 
300° C. No measurable weight loss due to absorbed water was found. In Fig. ||the accumulated 7-ray intensity from 
one of the '^^S activations is shown. The 7-line is well isolated on a low level of background counts. 

The following systematic uncertainties were combined by quadratic error propagation: Au standard cross section: 
1.5-3%, Ge-detector efficiency: 6%, 7-ray intensity per decay: 0.6% for the ^''^S and 0.1% for the ^^*Au decay, 
divergence of neutron beam: 2-12%, factor /;,: 1.5%, sample weight: <0.5%, and other systematic errors: 2%. 



IV. THEORETICAL ANALYSIS 



A. Reaction mechanism and models 



In nuclear reactions two extreme types of reaction mechanisms can exist: the compound-nucleus (CN) and the 
direct (DI) process. In the CN mechanism the projectile merges in the target nucleus and excites many degrees of 
freedom of the CN. The excitation proceeds by way of a multistep process and therefore has a reaction time typically 
of the order 10~^^s to 10"^" s. After this time the CN decays into various exit channels. The relative importance of 
the decay channels is determined by the branching ratios to the final states. In the DI process the projectile excites 
only a few degrees of freedom (e.g. single-particle or collective). The excitation proceeds in one single step and has 
a characteristic time scale of 10~^^ s to 10~^^ s. This corresponds to the time it takes the projectile to pass through 
the target nucleus; this time is much shorter than the reaction time of CN processes. 



B. Folding procedure 



The folding procedure is used for calculating the nucleon-nucleus potentials to describe the elastic scattering data 
and the bound states. This method has been applied successfully in describing many nucleon-nucleus systems. In the 
folding approach the nuclear density pA is folded with an energy and density dependent nucleon-nucleon interaction 

V,B 0, 111 



ViR) = XVf{R) - A / pA{r^v,s{E, pa, \R ~ n)dr 



(3) 



with R being the separation of the centers of mass of the two colliding nuclei. The normalization factor A is adjusted 
to reproduce the elastic scattering data the binding energies of the residual nuclei. The folding potentials of Eq. ^ 
were determined with the help of the code DFOLD |32|. 



C. Direct— capture model 



In thermonuclear scenarios the projectile energy is well below the Coulomb and/or centrifugal barrier. Consequently, 
the CN formation may be suppressed, because there are almost no CN levels that can be populated, especially in 
light, magic and far-off-stability nuclei. In this work the theoretical cross sections are calculated by only considering 
the DC-contributions. The theoretical cross section tr*'^ is obtained from the DC cross section a^'^ by ||3l[], Q 

a'^^Y.C"^^'^^'' ■ (4) 

i 

The sum extends over the ground state and excited states in the final nuclei, where the spectroscopic factors Si are 
known. The isospin Clebsch-Gordan coefficients are given by Ci . The DC cross sections af^ are essentially determined 
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by the overlap of the scattering wave function in the entrance channel, the bound-state wave function in the exit 
channel and the multipole transition-operator. The radial dependence of the wave functions in the DC-integral is in 
our case determined uniquely by the folding potentials. 



V. CALCULATIONS AND RESULTS 



In the folding approach the nuclear density pA for the stable nucleus ^^S was derived from the experimental charge 
distribution The normalization factor A for ^^S(n,7)^^S of the optical potential in the entrance channel was 

adjusted to fit the thermal (^^S+n)-scattering cross section of (1.1 ± 0.8) barn |Q. Even so this cross section is 
not determined well, we fitted our normalization factor to reproduce 1.1 barn. However, applying the same fitting 
procedure to the (^''S+n)-scattering cross section that is known much better ((1.52 ± 0.03) barn [^) we obtained 
almost the same volume integral per nucleon in the two cases (^''S+n: 501.7MeVfm'^; ^^S+n: 497.1 MeVfm'^). The 
imaginary part of the optical potential is small for the (^^S+n)-channel and can be neglected. For the exit channels 
the normalization constants A were adjusted to the energies of the ground and the excited states. The potentials 
obtained in this way ensure the correct behavior of the wave functions in the nuclear exterior. 

The spectroscopic factors for one-nucleon stripping of "^^S were determined from the most recent experimental 
■^^S(d,p)^^S-data ||3^. We also carried out shell-model calculations using a combined sd- and pf-shell with an effective 
nucleon-nucleon interaction derived by Warburton et al. jsoj. For these calculations the program OXBASH was 
used to calculate the wave functions and spectroscopic factors. The calculated spectroscopic factors are larger than 



those extracted from the experimental '^^S(d,p)'^^S-data (see Table III ) . This fact might be due to the strong truncation 
of the model space (only one neutron in the pf-shell for the negative parity states). Inclusion of more particle-hole 
excitations should give smaller spectroscopic factors. The masses and Q-values for the transitions to the different 



states of the residual nucleus ^^S were taken from experimental data ||3^, |39|. For the DC-calculations the code 
TEDCA @ was used. 

The level scheme of the relevant levels for '^^S is shown in Fig. ^. The cross section for the reaction '^^S(n,7)^^S 
obtained from the DC-calculation is compared with the experimental data from the thermal to the thermonuclear 
energy region in Fig. There are two types of El-transitions contributing for the transitions to the residual nucleus 
^^S. The first one comes from an s-wave in the entrance channel exciting the negative-parity states 3/2~ and l/2~ 
(see Table |l|). These transitions give the well-known 1/v-behavior (see Fig. ||). The second type of El-transition 
comes from an initial p-wave and excites the positive-parity state 3/2+ in the final nucleus. This transition has a 
v-behavior and can be neglected in the relevant energy range (see Fig. || and Table III). The El-transition to the 



7/2 -ground state of S can also be neglected, because of the higher centrifugal barrier necessary for the incoming 



d-wave (see Table [II ) . As can be seen from Fig. |j this contribution effects the deviation from an 1 / v-behavior of 
the cross section only above about 700 keV. 

The spin and parity assignments of the final states in ■^^S, the Q~values for the transitions to the different final 
states, the spectroscopic factors obtained from ■^^S(d,p)^^S and shell-model calculations are shown in Table III. Also 



in this table the calculated cross sections for •^^S(n,7)'^^S at 25.3meV, 25, 151, 176, and 218keVkeV using DC with 
the spectroscopic factors obtained from the (d,p)-data are compared with the experimental data. 

We have determined the thermonuclear-reaction-rate factor Na (crv) |^^. Since the cross section follows an l/w-law 
up to 150 keV we obtain a constant reaction-rate factor 

Na {av) = 2.56 x 10^ cm^ mole"^ s"^ . (5) 



VI. DISCUSSION 



Direct-capture calculations using the folding procedure can excellently reproduce the non-resonant experimental 
data for the capture cross section by the neutron-rich sulfur isotope ^^S in the thermal and thermonuclear energy 
region. The enhancement in the region of 176 keV comes from resonant contributions | |39t not considered in the 
DC-calculation. 

DC is also the dominant reaction mechanism for neutron capture by neutron-rich isotopes far-off stability occurring 
in the a-and r-process. For such isotopes the Q-value and therefore the excitation energy of the compound nucleus 
gets still lower, leading to a further substantial diminution of the level density of the compound nucleus. Thus, 
the DC-contribution becomes the dominating reaction mechanism. Nuclear-structure models are indispensable for 
extrapolating reaction rates to nuclei near and far-off the region of stability, because only a limited or no experimental 
information is available in this region. The DC cross reaction ''^S(n,7)^^S can be considered as a benchmark of different 
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nuclear-structure models (Shell Model, Relativistic Mean Field Theory, Quasi Particle Random Phase Approximation, 
Hartree Fock Bogoliubov Theory) for calculating neutron-capture cross sections by neutron-rich nuclei taking place 
in the a- or r-process . 

The s-process production of "^^S was recently discussed quantitatively by Schatz et al. but without a reliable 
'^^S(n,7) cross section. The s-process reaction network in the sulfur to calcium region contains (n,7), (n,p), and (n,Q!) 
reactions. The ■^^S production is mediated by the ■^^Cl(n,p)'^®S reaction from seed nuclei with mass numbers <36. But 
also seed nuclei >36 can contribute through the '^^Ar(n,Q;)'^^S reaction channel. Besides its formation the destruction 
of '^^S by the (n,7) reaction is important. A decrease in the '^^S(n,7)'^^S cross section leads to a corresponding increase 
of the abundance formed. As our measured ^^S(n,7)^^S value is by a factor of 1.8 smaller than the estimate of 
Woosley et al. the s-process abundance production of ^^S will be enhanced by this factor. The quantitative 
analysis requires also model parameters for the main and especially the weak s-process component. This information 
can be obtained from the analysis of the s-process beyond A=56 
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TABLE I. Sample characteristics and decay properties of the product nuclei S and Au 



Isotope 


Chemical 


Enriclimoiit 


Reaction 


Tl/2 


E-. 


Iiit(>risitv per dc;(;ay 






(%) 






(koV) 


(%) 




sulfur powder 


5.933 


^^'■S(n,7)'^^''S 


5.05 mill 


3103 


94.0±0.6 


i^^Au 


metallic 


100 


i«^Au(n,7)i'-'«Au 


2.69 d 


412 


95.50±0.096 
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TABLE II. Sample 


weights and experimental ^^S capture cross sections 






Mean neutron 


Mass of Au 


Mass of 


Mass of Au 


Irradiation 




O" 


Uncertainty 


energy 


front side 


sulfur 


back side 


time 




(^barn) 


statistical 


total 


(keV) 


(mg) 


(mg) 


(mg) 


(d) 






(%) 


(%) 


25 


16.80 


100.32 


16.84 


3.54 




192 


1.2 


11.3 




16.61 


95.32 


16.64 


1.70 




195 


1.1 


10.5 




16.57 


50.40 


16.51 


1.16 




169 


2.8 


8.3 




16.57 


50.40 


16.51 


1.70 




168 


2.3 


8.1 




15.71 


148.89 


15.73 


0.58 




191 


1.5 


13.9 




8.86 


50.11 


8.80 


2.70 




178 


1.0 


8.3 




17.15 


21.18 


17.22 


5.10 




190 


1.6 


7.5 










Average 




187±14 






151 ± 15 


16.38 


94.96 


16.71 


3.70 




85 


3.3 


9.5 




16.46 


150.19 


16.33 


6.90 




80 


1.8 


11.5 










Average 




81±7 






176 ± 20 


17.35 


94.94 


17.34 


2.40 




114 


7.2 


11.4 




17.35 


94.94 


17.34 


7.80 




128 


3.2 


9.4 










Average 




125±11 






218 ±23 


17.39 


150.00 


17.39 


3.20 




87 


6.3 


11.2 




17.40 


149.85 


17.38 


6.80 




74 


4.7 


11.7 










Average 




78±7 






TABLE III. Final states, Q-values, transitions and spectroscopic factors for the states of ^^S obtained from shell-model 


calculations and ^® 


S(d,p)^^S. Cross sections for 


^'^S(n,7)^'^S at 25.3 meV, 25keV, 151 keV, 176 keV, and 218 keV using DC with 


the experimental data. 
















Final Q-value 


Transition 


Spectroscopic factor 






Cross section 






state [MeV] 




shell 


(d,p)-reaction 


25.3 meV 25keV 


151 keV 


176keV 


218 keV 






model 


PI 


[mbarn] [/ibarn] 


[/ibarn] 


[^barn] 


[/ibarn] 


1" 4.303 


d ^ f 


0.91 


0.50 


0.0 


0.0 


0.1 


0.2 


0.2 


1" 3.657 


s ^ p 


0.86 


0.55 


157.1 


158.0 


64.3 


59.5 


53.5 


1+ 2.906 


p -> d 


0.07 


0.03 


0.0 


0.0 


0.0 


0.0 


0.0 


2.312 


s ^ p 


0.12 


0.03 


5.3 


5.3 


2.1 


2.0 


1.8 


1" 1.666 


s ^ p 


0.97 


0.47 


28.3 


28.5 


11.6 


10.8 


9.7 






Total cross section: DC 


190.7 


191.8 


78.1 


72.5 


65.2 




Total cross section: experiment 


150 ±30' 187 ±14 


81 ±7 


125 ± 11 


78 ±7 



Ref. li) 



FIG. 1. Scheme of experimental setup 

FIG. 2. Neutron spectra with mean neutron energies of 151(top), 176(middle), and 218 keV(bottom) at sample position. 
Using the proton energies and the metallic Li-targets (2.5 /im thickness) of the actual activation measurements they were 
calculated from TOF spectra measured under different angles with respect to the beam axis. 



FIG. 3. Accumulated intensity of the 3103 keV S 7-line from the activation with a 100.32 mg sulfur sample 



FIG. 4. Level scheme of ^^S 
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FIG. 5. Comparison of the DC cross section for ^^S(n,7)^'^S with the experimental data from thermal to thermonuclear 
projectile energies. The DC contributions for the different transitions to the final states of ^^S as well as the sum of all 
transitions (solid curve) is shown. The experimental data at the thermal energy have been taken from Ref . [35] . 
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Q=4.303 MeV 



36 



S+n 



E [MeV] 



2.638 



2.023 



1.992 
1.398 



1/2' 



7/2' 



3/2" 
3/2"^ 



0.646 



3/2' 



7/2 
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